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We perform photoluminescence experiments at 4 K on two different transition metal diselenide 
monolayers, namely MoSe 2 and WSe 2 in magnetic fields up to 9 T applied perpendicular to 
the sample plane. In MoSe 2 monolayers the valley polarization of the neutral and the charged 
exciton (trion) can be tuned by the magnetic field, independent of the excitation laser polarization. 

In the investigated WSe 2 monolayer sample the evolution of the trion valley polarization depends 
both on the applied magnetic field and the excitation laser helicity, while the neutral exciton valley 
polarization depends only on the latter. Remarkably we observe a reversal of the sign of the trion 
polarization between WSe 2 and MoSe 2 . For both systems we observe a clear Zeeman splitting for 
the neutral exciton and the trion of about ±2 meV at ^ 9 T. The extracted Lande-factors for 
both exciton complexes in both materials are g « —4. 


I. INTRODUCTION 


The two-dimensional (2D) semiconductor structures 
based on monolayer (ML) transition metal dichalco- 
genides (TMDCs) show very strong light-matter inter¬ 
action, with an absorption of the order of 10% per ML 
in the visible region of the optical spectrum [1—3]. This 
strong interaction can be explored in microcavity physics 
and for laser and light emitting diode (LED) devices. 
The optical properties of ML TMDGs are governed by 
strongly bound excitons (binding energy is of the order of 
0.5 eV) [4-12]. The absence of an inversion centre in the 
lattice of TMDG MLs together with the strong spin-orbit 
(SO) interaction in these materials leads to a coupling 
of carrier spin and fc-space valley dynamics. As a re¬ 
sult the circular polarization {a'^ or a~) of the absorbed 
or emitted photon can be directly associated with selec¬ 
tive exciton generation in one of the two non-equivalent 
AT-valleys: or K~, respectively [13-19]. Excitonic 

resonances do not only dominate single- and two-photon 
absorption, but also strongly influence second harmonic 
generation (SHG) [11] and interactions with plasmons 
[20]. The exact symmetry, degeneracy and energy spac¬ 
ing of the exciton levels is still under debate. To address 
the spin- and valley-dependent hne structure of excitonic 
levels we perform magneto-optical experiments on ML 
TMDGs, where the magnetic held lifts the underlying 
degeneracies and extract the effective Lande factors, a 
technique successfully used in the past to study the band- 
structure and excitonic effects of semiconductors [21]. 

Very recently lifting of the AT-valley degeneracy 
through the application of a magnetic held perpendic¬ 
ular to the monolayer plane has been reported for MoSe 2 
[22, 23] and WSe 2 [24, 25]. An energy splitting Az 
between the cr^ and cr” polarized photoluminescence 
(PL) components has been reported, which correspond 
to inter-band recombination in the K'^- and A'“-valley, 
respectively [13-17, 26]. The splitting Az increases 
roughly linearly with the applied magnetic held, although 
the values reported for ML WSe 2 differ by a factor of up 


to 3 [24, 25]. The aim of this comparative study is to 
perform experiments on MoSe 2 and WSe 2 monolayers in 
magnetic helds up to [S^] = 9 Tesla in the same set-up, 
and to measure the neutral exciton and trion emission 
energy and polarization. We vary the laser excitation 
energy and helicity, which allows us to distinguish be¬ 
tween the laser induced and magnetic held induced valley 
polarization, which show very different dependencies on 
the applied magnetic held when comparing MoSe 2 with 
WSe 2 monolayers. 


II. EXPERIMENTAL SET-UP AND SAMPLES 

WSe 2 and MoSe 2 hakes are obtained by exfoliation [27] 
of a bulk crystal on a Si02/Si substrate. The ML regions 
are identihed by optical contrast and very clearly in PL 
spectroscopy. Experiments at T = 4 K and in magnetic 
helds up to ±9 T have been carried out in an ultra-stable 
confocal microscope developed for spectroscopy on single 
semiconductor quantum dots (typical diameter of 20 nm) 
[28, 29]. We can infer that the mechanical movement of 
the detection and excitation spot on the sample due to 
the strong applied magnetic helds is at most in the tens 
of nm range in our set-up. The detection spot diameter 
is about 700 nm. The sample is excited either by a He- 
Ne laser (1.96 eV) or by a tunable continuous wave Ti:Sa 
laser. The average laser power is in the /iW range, in the 
linear absorption regime. The PL emission is dispersed 
in a double-monochromator and detected with a Si-GGD 
camera. The spectral resolution of this detection system 
is « 20 /ieV. The circular PL polarization Pc is dehned 
as 

Pc = (/.+ -/.-)/(/.++/.-), (1) 

where I„+{Icr-) denotes the intensity of the right (cr+) 
and left ((T“) circularly polarized emission. Light is a'^ 
(right circularly polarized) if electric held vector rotates 
with time clock-wise provided one looks along the light 
propagation axis. Similarly the linear PL polarization 
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FIG. 1: Data for monolayer MoSe2; excitation energy i^Laser = 1-96 eV (a) PL spectra at Bz = -1-8 T for (black) and 
(j~ (red) polarized detection using a'^ polarized laser excitation, (b) same as (a) but at Bz — 0 T. (c) same as (a) but at 
Bz = —8 T. (d) Splitting between the and a~ polarized PL components for the trion (black squares) and the X°(blue circles) 
as a function of magnetic field, the function —A^bBz is shown for comparison (red line), (e) Schematics of the experimental 
geometry, (f) Polarization of the PL emission of the X'^ as a function of Bz- (g) Polarization Pc of the PL emission of the trion 
as a function of Bz using cr^ (black squares) and a~ excitation (red squares). 


in the fixed {xy) axes lying in the sample plane writes 
Pi = {Ix - ly)/{Ix + ly) with Ix{Iy) the X and y linearly 
polarized emission components. 


III. MAGNETO-OPTICS IN MONOLAYER 
MoSe2 

First we discuss the experimental results for MoSe 2 
MLs. At zero magnetic field, we observe in the PL spec¬ 
trum two sharp emission features with a full width at 
half maximum (FWHM) of 10 meV. In accordance with 
previous reports [6], the low energy emission at 1.63 eV 
is attributed to the charged exciton (trion) recombina¬ 
tion and at 1.67 eV we record the neutral exciton X° 
emission, see Fig. lb. In our measurements the circular 
polarization degree Pc of the trion and X° emission de¬ 
pends only very little on the excitation laser polarization 
(cr+ or a~) for the excitation laser energy of 1.96 eV used 
in Fig. lb. The maximum Pc generated was of the order 
of 5% for the excitation energy range investigated [30], 
much lower than the PL polarization achieved for circu¬ 
lar laser excitation in M 0 S 2 [14-17, 26] and WSe 2 MLs 
[11, 18, 31]. Below we discuss the magnetic field depen¬ 
dence of the PL emission, whose energy and polarization 
are essentially independent of the excitation laser polar¬ 
ization, in stark contrast to the results obtained in ML 
WSe 2 discussed in Sec. IV. 

Next we discuss the changes observed in the PL emis¬ 


sion when applying a magnetic field B^, perpendicular to 
the layer plane, i.e. along the z-direction, that is also the 
light propagation axis (Faraday geometry). Comparing 
Fig. la at = 4-8 T and Fig. lb at = 0 we ob¬ 
serve two main differences: (i) For Bz = 4-8 T the a'^ 
polarized emission is more intense than the a~ polarized 
component, for both the trion and X°; (ii) The a~^ po¬ 
larized emission is shifted to lower energy compared to 
the a~ polarized component. When applying a field of 
Bz = —8 T in Fig. Ic, these results are reversed, i.e. a~ 
emission becomes more intense and lower in energy com¬ 
pared to a~*' in agreement with time-reversal symmetry. 
In Fig. Id we plot the full magnetic field dependence of 
the energy splitting 

Az = = gysBz, ( 2 ) 

from Bz = —9 T to 4-9 T with g being the effective g- 
factor and fiB being the Bohr magneton. The Zeeman 
splitting is extracted by fitting the trion and X*^ emis¬ 
sion spectra with Lorentzians. We observe a clear linear 
dependence both for the trion and X° splitting on the ap¬ 
plied field. This indicates that eventual diamagnetic and 
higher order contributions are identical for the upper and 
lower Zeeman branch within our experimental resolution. 
The slope is AzjBz = —220±10 /reV/T for the X° which 
corresponds to an exciton g-factor of gx^ = —3.8 ± 0.2. 
For the trion the slope is —226 ± 10 /ieV/T which cor¬ 
responds to a (jf-factor of gx = —3.9 ± 0.2. This corre¬ 
sponds to a maximum Az = —2 meV at 9 Tesla. The 
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main experimental uncertainty for Az comes from even¬ 
tual changes of the overall shape of the PL emission due 
to imperfections in the optical set-up and sample inho¬ 
mogeneities, as the recorded shifts are smaller than the 
linewidth. The Lande-factors for trions and the X° in 
ML MoSe 2 extracted from our date are close to the val¬ 
ues reported in [22, 23]. 

It is worth mentioning that trions that form with the 
excess electron in the same or different valley, with re¬ 
spect to the photo-generated electron hole pair, are sep¬ 
arated in energy due to the strong Coulomb effects and 
the zero-field splitting of conduction band states in each 
valley, cf. [32]. However, we do not observe any fine struc¬ 
ture splitting for the trion emission in the investigated 
sample. Hence, the measured gx represents the global 
magnetic field induced energy shift of the trion emission. 

We now discuss the circular polarization of the ob¬ 
served emission summarized in Figs. If and Ig. For the 
trion. Pc increases from zero to 30% when the magnetic 
field is increased from zero up to 9 T. This strong po¬ 
larization is reversed when the direction of the applied 
magnetic field is reversed. Optical valley initialization at 
higher fields is not at the origin of this effect: For cr+ 
polarized laser excitation we observed exactly the same 
polarization increase with field as for cr“ polarized laser 
excitation. We conclude that the PL polarization is the 
result of magneto-induced spin/valley polarization build¬ 
up during the PL emission time, which has been deter¬ 
mined to be in the ps-range [30]. Observing an increase 
in the emission polarization in applied magnetic fields is 
a very common observation in semiconductors [33]. For 
ML MoSe 2 , where excitons have a strong binding energy 
[8] , this observation is very surprising taking into account 
the relevant energy scales. At 9 Tesla, we have induced 
a splitting between the valley Zeeman levels of 2 meV, 
eventually comparable to the spin splitting in the conduc¬ 
tion band, predicted to be in the meV range [34, 35]. The 
polarization of the neutral exciton follows a similar trend 
as the trion. For the neutral exciton, the applied field 
Bz can dominate the long-range electron-hole Coulomb 
exchange interaction, and hence suppress valley depolar¬ 
ization when has a larger amplitude than the effective 
transverse field associated to exchange effects [32, 36-38]. 


IV. MAGNETO-OPTICS IN MONOLAYER WSez 

In order to draw more general conclusions on valley 
properties in applied magnetic fields, we compare our re¬ 
sults in ML MoSe 2 with a very well characterized system, 
namely ML WSe 2 [11, 18, 31] investigated with the same 
experimental set-up. Therefore all experimental uncer¬ 
tainties are comparable. In Fig. 2b the PL emission of 
ML WSe 2 is plotted, and just as in the case of ML MoSe 2 , 
we observe two well defined emission lines with a typi¬ 
cal FWHM of 10 meV. The higher energy transition at 
1.75 eV is identified as the neutral exciton X° recom¬ 
bination, as valley coherence can be generated [18], see 


TABLE I: comparison of Lande factors 



9x0 

X° slope 

gx 

Trion slope 



[peV/T] 


[geV/T] 

ML MoSe 2 

-3.8 ± 0.2 

-220 ± 10 

-3.9 ±0.2 

-226 ± 10 

ML WSe2 

-3.7 ±0.2 

-214 ± 10 

-4.4 ±0.2 

-254 ± 10 


Fig. 2e. The lower energy transition at 1.72 eV stems 
from the radiative recombination of the charged exciton 
(trion). The and cr~ polarized PL components corre¬ 
spond to carrier recombination in the K'^ and K~ valley, 
respectively. At zero magnetic field the emission spectra 
in a~^ and a~ polarizations are exactly the same due 
to the time-reversal symmetry and corresponding degen¬ 
eracy of the valley levels. In a magnetic field of 9 T 
this valley degeneracy is lifted, and the energy difference 
Az = Eq. (2), is typically —2 meV, see 

Fig. 2a, very similar to the results on MoSe 2 presented 
in Fig. la. In a field of = —9 T, the (t+ polarized 
component is now at higher energy, corresponding to a 
positive Zeeman splitting of +2 meV. We plot in Fig. 2d 
the full magnetic field dependence of the Zeeman splitting 
for the trion and X°. The experimental points are very 
close to a straight line, which allows us to extract for the 
X° a slope oi Az/B^ = —214±10 ^eV/T, corresponding 
to the neutral exciton ^-factor gxo = —3.7 ± 0.2. For 
the trion the slope is —254 ± 10 fieV/T, resulting in its 
(/-factor gx = —4.4 ± 0.2. 

The similarities between the measured Lande-factors 
in MoSe 2 and WSe 2 are striking, both with respect to 
their signs and amplitudes, see table 1. Comparing with 
the very recent literature data on ML WSe 2 , the g- 
factors for the neutral exciton extracted for our sample 
are smaller than reported in reference [25] and larger than 
in reference [24]. Note, that as for MoSe 2 case we do not 
resolve any fine-structure of the trion emission in WSe 2 . 

Although the Zeeman energy evolution of the both ML 
materials MoSe 2 and WSe 2 shows close similarities, we 
will see below that the evolution of the valley polarization 
as a function of the applied magnetic field is completely 
different. The polarization measurements on WSe 2 are 
all performed with non-resonant HeNe laser excitation. 
First, we plot the evolution of the X° PL polarization 
under polarized pumping, that corresponds to the val¬ 
ley polarization generation via the chiral optical selec¬ 
tion rules [13, 14], as a function of the applied magnetic 
field. It can be seen in Fig. 2f that the X° polarization 
is essentially independent of the applied field. For a lin¬ 
early polarized excitation laser, we find Pc ^ 0 for all 
applied field values. Changing to a~^ circularly polarized 
excitation, we record at zero Tesla Pc « 40%. This value 
remains practically constant in applied fields from —9 T 
to ±9 T. This indicates that over the entire magnetic field 
range, the X° emission polarization is determined by the 
initially, optically created valley polarization rather than 
induced by the magnetic field. During the X° PL emis¬ 
sion time no valley/spin relaxation occurs, which might 
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FIG. 2: Data for monolayer WSe2; excitation energy i?Laser = 1-96 eV (a) PL spectra at = -|-9 T for (black) 
and a~ (red) polarized detection using polarized laser excitation, (b) same as (a) but aX Bz = Q T. (c) same as (a) but 
at Bz = —9 T. (d) Splitting between the cr^ and a~ polarized PL components for the trion (black squares) and the X®(blue 
circles) as a function of magnetic field, the function —AfisBz is shown for comparison (red line), (e) Linear polarization of the 
X° emission as a function of Bz- (f) Polarization of the PL emission of the X° as a function of Bz for three different laser 
polarizations, (g) Polarization of the PL emission of the trion as a function of Bz for three different laser polarizations. 


simply be a consequence of the extremely short PL emis¬ 
sion time in the few ps range measured for this material 
[31]. Exciting preferentially the K~ valley with a cr“ po¬ 
larized laser, results in Pc « —40%, again independent of 
the applied magnetic field. This behaviour for the X*’ in 
ML WSe 2 is in stark contrast to the observations in ML 
MoSe 2 . For the former, the optical valley initialization 
determines the PL polarization, for the latter the applied 
magnetic field direction and amplitude allow to control 
the valley polarization. Comparing with the recent liter¬ 
ature, the X° polarization in the ML WSe 2 sample inves¬ 
tigated by Aivazian et al. [24] was slightly more sensitive 
to the applied magnetic field. 

Contrary to the X° in our ML WSe 2 sample, the trion 
polarization in this ML can be controlled via the applied 
magnetic field, shown in Fig. 2g. For linearly polarized 
laser excitation, the trion PL polarization for Bz = 0 is 
absent. Pc = 0, but increases at —9 T to -|-30%, and at 
-1-9 T to —30%. This can be directly compared to the 
trion in ML MoSe 2 , that also started off at zero field 
with zero polarization. When comparing the trion polar¬ 
ization evolution in ML MoSe 2 (Fig. Ig) and ML WSe 2 
(Fig. 2g) we note very contrasting behaviour: Applica¬ 
tion of a positive magnetic field, results in a strong, pos¬ 
itive Pc in MoSe 2 in contrast to the strong, but negative 
polarization created in WSe 2 . Using cr+ polarized exci¬ 
tation for the trion, results in high circular polarization 
Pc = 50% already at Bz = 0. At a magnetic field of 
Bz = —9 T this polarization is increased to Pc = 68%, 


at Bz = -1-9 T we find Pc = 34%. For the trion using 
a~^ excitation, we find a similar trend as under linearly 
polarized excitation: the PL polarization decreases for 
positive Bz and increases for negative Bz- The trion 
polarization depends on both the excitation laser polar¬ 
ization and the applied magnetic field. This latter de¬ 
pendence might be linked to the fact that the PL emis¬ 
sion time of the trion in ML WSe 2 is longer than for 
the X° [31], allowing for polarization relaxation to oc¬ 
cur before radiative recombination takes place. As can 
be seen in Fig. 2g, the magnetic field induced change in 
PL polarization is also observed for a~ excitation, here 
we find at Bz = —9 T Pc = —31%, while for zero field 
we record Pc = —52% (as expected when switching from 
a~^ to a~ excitation) evolving towards Pc = —69% for 
Bz = —9 T. The three different measurement series pre¬ 
sented in Fig. 2g each confirm that the PL polarization 
of the trion in ML WSe 2 strongly decreases with the ap¬ 
plied magnetic field, whereas, in stark contrast, the trion 
PL polarization in ML MoSe 2 strongly increases as a 
function of the applied Bz- 


V. DISCUSSION 

The emission energies and polarizations are deter¬ 
mined by several factors. Most importantly, the Zeeman 
splittings of exciton and trion in the emission spectra 
are governed by an interplay of spin splittings of conduc- 
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tion and valence band states, while the polarization is 
governed by (i) selection rules at optical transitions, (ii) 
occupancies of the spin/valley states and possible spin 
relaxation processes. Here we briefly discuss theoretical 
approaches to evaluate Zeeman splittings and polariza¬ 
tion of emitted radiation in the context of experimental 
data presented above. 


are the interband matrix elements of the coordinate op¬ 
erator. In derivation of Eq. (6) we made use of the com¬ 
pleteness relation for the Bloch amplitudes and took into 
account that = 0. Taking into account that 

E^)mo\ with E„, Em being the en¬ 
ergies of corresponding bands, Eq. (6) can be rewritten 
as [40] 


A. Zeeman effect in two-dimensional crystals 

We recall that the electron wavefunction in the two- 
dimensional crystal can be recast, in accordance with the 
Bloch theorem, as 


i^nq{r) 


^iqr 




( 3 ) 


where n enumerates bands (including electron spin state) 
and q is the quasi-wavevector, Unq{r) is the periodic 
amplitude normalized per volume of the unit cell, Hq: 
Ifio \'^nq\'^d'f = ^^ 0 ) and S is the macroscopic normaliza¬ 
tion area. 

The bare electron Zeeman effect is described by 


■Ho = • s, (4) 

where pB = \e\h/{2mQc) is the Bohr magneton, mo and 
go = 2 are the free electron mass and Lande factor, and 
s is the spin operator (t/ 2, cr being the vector composed 
of the Pauli matrices. In crystals, in addition to (4) the 
orbital contribution to the Zeeman splitting should be 
taken into account. It is related to the orbital momentum 
of the electron [39] , 


Hi = pbB ■ L, (5) 

where L = hr^[r x p] is the angular momentum opera¬ 
tor, and p = —i/iV is the electron momentum operator. 
The contribution (5) is known to be important for well 
characterized semiconductors such as GaAs [21, 40]. 

We recall that in MoSe 2 and WSe 2 MLs the direct band 
gaps are realized at the edges of the Brillouin zone char¬ 
acterized by the wavevectors . The bands aX q — 
are non-degenerate and can be characterized by a certain 
spin projection Sz = ±1/2 onto the sample normal. The 
time reversal symmetry couples = ±1/2 states in the 
valley with Sz = ±1/2 states in the K~ valley. In 
what follows we will be interested in the Zeeman effect in 
magnetic field B \\ z. Correspondingly, it is instructive 
to present the diagonal matrix element of operator at 
q = K+ or K- as 

{tpn\^Lz\i>n)± = 

m^n 

Here (a = cc, y) are the matrix elements of the elec¬ 
tron momentum and = i^o ^ Iqi , K (du m /dqa)dr 


^ y plm{K^)P^mniK^) " plm{K^)P%.niK^) ^ 
•mo ^ Em- En 

m^n 

( 7 ) 

Equations (4), (5) and (7) can be used to evaluate Zee- 
man splittings both in the k p and tight-binding models, 
as detailed in the next two subsections. 


B. kp-theory 

The multiband k p model was formulated for TMDCs 
in Refs. [41-43]. The effective Hamiltonians describing 
the states in the vicinity of edges of the Brillouin 
zone have the form 


/Ec+2 

76 fc- 

74 /C+ 

0 

l&k+ 

Ec 

73 fc- 

l3k+ 

14 k- 

l3k+ 

Ey 

12 k- 

\ 0 

75 fc- 

l2k+ 

Ey-3 

/Bc+2 

lek+ 

14 k- 

0 

lek- 

Ee 

l3k+ 

75 fc- 

, I4k+ 

73 fc- 

Ey 

l2k+ 

\ 0 

l5k+ 

to 

1 

Ey-3 


(8a) 


(8b) 


Here k± = ± iky are the cyclic components of the 

electron wavevector reckoned from the points, k± = 
— q\ the parameters 73 ... 75 are related to the inter¬ 
band momentum matrix elements, and certain conven¬ 
tion about the phases of the Bloch functions is assumed. 
The symbols c ± 2 and v — 3 denote excited conduction 
and deep valence bands. Such a model was shown to al¬ 
low for adaquate description of the spin-orbit coupling 
and trigonal symmetry effects in TMDCs [41-43]. In the 
general theory besides fc-linear off-diagonal terms, effec¬ 
tive k ■ p Hamiltonian includes diagonal quadratic in k 
contributions resulting (i) from the second-order k ■ p 
coupling with distant bands and (ii) from bare electron 
dispersion h'^k^/2mo [44]. Usually these contributions 
are comparable and should be included simultaneously. 
Diagonalizing Hamiltonians (8) in the second order in off- 
diagonal terms we obtain for the electron effective masses 
of the conduction and valence bands, respectively, 

1 1111 

— “ —7 ^- 7 ^-'- 7 ’ 

TOc m* m/ mo m/ 

1 1111 

- “-7 "I- 7 "I-'-7 

my m* m'y mo m" 


(9a) 

(9b) 
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Here subscripts c, v denote the corresponding bands and 
we use the electron representation, m* = — 

Ey)/{2j^) is the two-band effective mass, the terms 

2 2 

li + ^6 

Ec — Ey -3 Ec — Ec^2 




± = 1 ( ^2 ll 

m'y \Ey - Ey _3 Ey-Ec+2 

result from the mixing described by the Hamiltonians (8) 
and 1/™" contain above mentioned contributions 

from remote bands [not included in Eqs. (8)]. 

According to the general theory [44] the magnetic field 
within the k ■ p scheme is included (i) by adding the 
bare Zeeman effect in the form of Eq. (4) and (ii) by 
replacing k in Eqs. (8) by fc — {ejch)A, where e = —|e| 
is the electron charge, A is the vector potential of the 
magnetic field. The calculation in the first order in 
yields the effective g-factors of electrons in valleys: 



9c,V = 9c,V = 2 5 ” , 


The term 2 in Eq. (11) arises from the bare Zeeman effect, 
Eq. (4) , while g°'^y result from the k ■ p-mixing 


orb ^ 4 too f ll _ ll 

^2 \Ec-Ey Ec-Ey_3 


7i 


Ec — E, 


c+2 


(12a) 


orb ^ f ll ll ll \ 

\ Ey-Ey Ey-Ey_3 

(12b) 

Equations (12) can be derived from Eqs. (5) and (7), 
taking into account only bands c, c -I- 2, u and w — 3. 
Additional contribution to g°^y , namely, Ag°‘^y may arise 
allowing for other distant bands in Eq. (7). 


C. Tight-binding theory 

In the tight-binding approximation the Bloch function 
in Eq. (3) is presented as a linear combination of atomic 
orbitals (j)'j{r) (generally orthogonalized [45]) in the form 

^n,{r) = ^ - Rc,i). (13) 

a,1,3 

Here a enumerates types of atoms (a = Metal or Chalco- 
gen), I enumerates atoms of a given type, j runs through 
the set of orbitals taken into account at a given atom, 
Ra I are the positions of the atoms in the two-dimensional 
lattice and (n, q) are the coefficients. The 

tight-binding Hamiltonian acting in the space of coeffi¬ 
cients Cji contains diagonal energies, i.e. energies of or¬ 
bitals, and hopping matrix elements mixing coefficients 
C^j and Cp j, for different atoms {I ^ I'). There are sev¬ 
eral tight-binding models developed for TMDC monolay¬ 
ers [46]: (i) the 3-band model [35] which takes into ac¬ 
count only three d-orbitals of Metal, namely, 


and dxy, and includes up to three nearest neighbours 
to reproduce density functional theory (DFT) calcula¬ 
tions of electron/hole dispersion over the whole Brillouin 
zone, (ii) the 7-band model [47] which includes three d- 
orbitals of Metal, two p-orbitals for each Chalcogen and 
only nearest neighbour interactions (this model involves 
as well overlap matrix elements since it is developed for 
non-orthogonalized orbitals); (iii) the 11-band model of 
Ref. [48] which accounts for five d-orbitals of Metal and 
three p-orbitals of each Chalcogen and applies the near¬ 
est neighbour approximation; (iv) the 27-band model of 
Ref. [49] which takes into account sp^d^ states for each 
atom and provides rather high accuracy. 

The diagonalization of the tight-binding Hamiltonian 
provides the energy dispersion in the whole Brillouin 
zone. In the vicinity of points the tight-binding 
Hamiltonian can be expanded in power series in the elec¬ 
tron wavevector (referred to the Brillouin zone edge) and 
can be partially diagonalized. The resulting effective 
Hamiltonian has a form similar to the k ■ p Hamilto¬ 
nian (with different number of bands depending on the 
tight-binding model) but contains diagonal contribu¬ 
tions. It is convenient to interpret the latter as an effec¬ 
tive k ■ p contribution of “the bare electron mass and re¬ 
mote bands”. The off-diagonal elements (linear in k) con¬ 
tain coefficients which can be interpreted as inter-center 
contributions to the momentum matrix elements [50-53] . 
For example, the 3-band model of Ref. [35] yields a 3 x 3 
effective Hamiltonian which includes the c, v and c -I- 2 
bands in notations of the k ■ p Hamiltonian (8). The 7- 
band model of Ref. [47] yields an additional conduction 
band with the same symmetry as the c-band. 

Within the tight-binding approach the magnetic field 
is also included in a two-fold way. First, the phase of 
the hopping matrix elements is modified by including the 
vector potential as follows 


ie 

ch 


A{r)dr. 


(14) 


This procedure is equivalent to the k ^ k — {e/ch)A 
replacement in the k ■ p Hamiltonian. Second, the intra- 
atomic contribution should be included. It contains the 
spin part, Eq. (4), as well as possible orbital contribu¬ 
tion. The latter should be carefully calculated using the 
restricted basis of atomic orbitals, the completeness rela¬ 
tion for this restricted basis, and an analogue of Eq. (6) 
in the form 


(15) 

i 

Note, that the summation in Eq. (15) extends over the or¬ 
bitals used in Eq. (13) rather than over all possible states 
of an isolated atom, otherwise the problem becomes ill- 
defined. 

In the particular cases of 3-, 7- and 11-band tight- 
binding models, the intra-site values of the momen¬ 
tum operator vanish due to symmetry reasons because 














7 


for any given atom either only p-shell or d-shell states 
are included into the model. Therefore, in Eq. (15) 
(^“(r)|liL^|^“(r)) = 0 and the intra-site contribution to 
the Zeeman effect vanishes. In this case we obtain within 
the tight-binding method Eq. (7) for the orbital mo¬ 
mentum and recover expressions analogous to Eq. (11), 
Eqs. (12) from k ■ p theory with no additional intra¬ 
center contributions in contrast to Refs. [22-24], see also 
Ref. [54] . Additional intra-site contributions may arise in 
the advanced tight-binding formalism of Ref. [49] where 
orbitals of different symmetry are included for a given 
atom, in which case the momentum operator matrix ele¬ 
ments between these orbitals may become non-zero. 


D. Zeeman splittings of direct excitons and trions 


The neutral exciton radiative decay involving emission 
of a'^ or a~ photons results from the recombination of a 
Coulomb-correlated electron-hole pair in the or K~ 
valley, respectively. The Zeeman splitting of X° is, in 
accordance with Eq. (2), given by 

^ - 9v^ - (sf ^ - 9^^)] 9bB^ . (16) 

Note that hereafter we neglect the renormalization of 
^-factor due to the Coulomb effects and band non- 
parabolicity. Making use of Eqs. (11), (12) the bright 
exciton (^-factor is given by 


9x — 9c 9v — 2 



mo \ 

m'y) 


+ Ag”*’ 


(17) 


where terms in parenthesis are calculated within the 
framework of Hamiltonian (4) and 4-band Hamito- 
nian ( 8 ) , and the contribution results form 

the remote bands not accounted for by Eq. ( 8 ). At this 
stage precise measurements of the conduction band effec¬ 
tive masses have not been reported yet. First measured 
values for the valence band effective masses can be ob¬ 
tained from angle-resolved photoemission spectroscopy 
(ARPES) [3, 55]. Equation (17) can also be represented 
in the form 


_ 2 (— + — - — - —^ 

\ TTic rriv to " to " J 

+ Ag^/'^ - AC" , ( 18 ) 


where the contributions to the Lande factor and effec¬ 
tive masses resulting from remote bands are explicitly 
present. Note that in the two-band approximation the 
bright exciton g-factor exactly vanishes because in this 
approximation the conduction- and valence-band elec¬ 
tron effective masses are given by toq /to^^ = 1 - 1 - toq /to* , 
mo/rriy^ = 1 — mojm* and, therefore. 


= 4-2 


rriQ 


mp 

2b 


TO 


= 4-2x2 = 0. 


The term 4 arises from inclusion of I/toq in the inverse 
effective masses for the conduction and valence band elec¬ 
trons. It should be emphasized that the inclusion of free 
electron dispersion I/toq in Eqs. (9) and neglecting the 
terms 1 /to", 1 /to/, Ag°‘^^, and Ag°'~^ contributed by re¬ 
mote bands cannot formally be justified within the k ■ p 
method. Those contributions should be estimated from 
experimental data or evaluated via more advanced k ■ p 
schemes or atomistic approaches. 

Since in MoSe 2 and WSe 2 the conduction- and valence- 
band effective masses are close in absolute values but 
have opposite signs, in accordance with Eqs. (17), (18), 
the orbital contribution to the exciton ^-factor stems 
from remote bands. We note that neither Eq. (17) with 
Ag”'^ — Ag”'^ = 0 nor Eq. (18) with to/, to/ —>■ oo can 
satisfactory describe the experiment in respect of both 
the sign and the magnitude of g-factor. Hence we con¬ 
clude that the remote band contributions are important. 
For instance, if we add into consideration one more (dis¬ 
tant) band with the same symmetry as the conduction 
band with band edge energy A/ and the matrix element 
oi k ■ p interaction with the valence band 73 we obtain 
an additional contribution to the valence band 5 -factor 
Agy = Amp^'^/[fi^{E'^ — Ey)], while Ag^. = 0 for the sym¬ 
metry reasons. Further experimental studies and theo¬ 
retical modelling are therefore needed to elucidate the 
values of Zeeman splitting. 

The Zeeman splitting for the bright trion is also given 
by Eq. (16). This is because the optical recombination 
involves charge carriers of opposite signs in the same 
valley, while the spin/valley state of the third carrier 
is not changed. The difference of measured values for 
excitons and trions can be attributed to the Coulomb- 
induced renormalization of 5 -factors due to band non- 
parabolicity. 


E. Polarization of emission 

The polarization of the PL emission is governed by the 
selection rules and occupancies of spin/valley states of 
the carriers and Coulomb complexes. In case of MoSe 2 , 
neither optical orientation (circular polarization of photo¬ 
luminescence for circularly polarized excitation at B = 0) 
nor exciton alignment or valley coherence (linear po¬ 
larization at linearly polarized excitation) are observed 
[30]. This allows us to assume that in this material 
the spin/valley relaxation of both neutral excitons and 
trions (or individual carriers) is fast compared to the 
PL emission time. In the presence of a magnetic field 
the magneto-induced circular polarization seems to result 
from the preferential occupation of the lowest Zeeman 
state of the exciton or trion. It is indeed consistent with 
experiment where, for > 0 , the state emitting the a'^- 
polarized photons has the lower energy and is dominant 
in the photoluminescence so that Pc{Bz > 0 ) > 0 , see 
Fig. 1. Note, however, that the full thermalization to the 
lattice temperature does not occur because at IB^] « 9 T 


the Zeeman splitting exceeds by far the temperature ex¬ 
pressed in the energy units, while experiment demon¬ 
strates only I Pc I 30 %. The effective spin temperature 
deduced from the experiment is about Tspin = 30 K. 

The situation is more complex for WSe 2 where the ex¬ 
periment shows a substantial optical orientation, Fig. 2f, 
and neutral-exciton alignment. Fig. 2e, even at = 0. 
In this case the spin/valley relaxation time is compa¬ 
rable to the lifetime of excitations and, in an applied 
magnetic field, thermalization may not occur. Addition¬ 
ally, as compared to MoSe 2 , WSe 2 is characterized by 
the opposite sign of zero-field spin splitting of the con¬ 
duction band [34, 35]. Hence, for the exciton and trion 
ground states the direct intra-valley optical transitions 
are spin-forbidden. These states can be manifested in op¬ 
tical spectra due to indirect (e.g. phonon-assisted) tran¬ 
sitions similarly to the case of Carbon nanotubes [56]. 
It follows from the symmetry considerations that in this 
case the polarization is reversed [57]. Further exper¬ 
imental data using gated devices where positively and 
negatively charged excitons can be clearly distinguished 
and theoretical analysis are needed to clarify this issue. 

Finally we emphasize that the above analysis is based 
on the perturbative treatment of the Coulomb interac¬ 


tion while, in TMDC MLs, the exciton binding energy 
amounts to > 0.5 eV and, hence, is comparable to the 
band gap. Theoretical estimates show that depending 
on the parameters of the materials and dielectric envi¬ 
ronment the binding energy is so large that the Is state 
emission could be in the infra-red range and, moreover, 
the exciton ground state could even collapse [58, 59]. In 
the former case the optical transition 2p —)■ Is could be 
relevant as well, its Zeeman splitting and polarization de¬ 
serve further study. In the latter case the ground state of 
the system could be strongly renormalized and an exci- 
tonic insulator could be formed [60-63]. Its polarization 
and magnetic field properties should also be studied in 
future works. 
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